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Abstract 

First results are presented from an analysis of data from the DRIFT-IIa and DRIFT- 
Ilb directional dark matter detectors at Boulby Mine in which alpha particle tracks 
were reconstructed and used to characterise detector performance — an important 
step towards optimising directional technology. The drift velocity in DRIFT-IIa was 
59.3 lb 0.2 (stat) ±7.5 (sys) ms~^ based on an analysis of naturally-occurring alpha- 
emitting background. The drift velocity in DRIFT-IIb was 57 it 1 (stat) ± 3 (sys) 
ms~^ determined by the analysis of alpha particle tracks from a ^^''Po source. 3D 
range reconstruction and energy spectra were used to identify alpha particles from 
the decay of '^'^'^Rn, ^iSpo, ^aoj^^ and ^i^Po. This study found that (22 ± 2)% of 
^""^^Po progeny (from ^^^Rn decay) are produced with no net charge in 40 Torr CS2. 
For ^^^Po progeny (from ^^'^Rn decay) the uncharged fraction is (100^35)%. 
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1 Introduction 



The DRIFT project aims to develop and operate tlie first underground time 
projection chamber (TPC) array suitable for observing and reconstructing 
WIMP-induced nuclear recoil tracks with enough precision to provide a signa- 
ture of the local WIMP Galactic halo [1,2,3,4]. The development of the DRIFT 
Negative Ion TPC (NITPC) detector concept is covered in [5,6,7,8,9,10,11,12]. 
A full desription of the latest detector design, DRIFT-II, including installa- 
tion and operation of the first modules, DRIFT-IIa and DRIFT-IIb, at Boulby 
Mine is covered in [12]. The purpose of this paper is to bring together results 
concerning the behaviour of DRIFT-IIa and DRIFT-IIb as evidenced by ob- 
servation of various classes of alpha particle events. The results in this paper 
are complementary to those covered in [13], where the alpha decay of radon 
progeny were found to produce recoil-like events. 

Although alpha particle events are, ultimately, an unwanted background in a 
rare event detector, their presence in DRIFT allows investigations of event dis- 
crimination, track reconstruction, directional sensitivity and other data anal- 
ysis capabilities. Determination of the accuracy of ion transport codes (e.g. 
SRIM2003 [15]) is also possibleQ. Alpha particles of a few MeV in 40 Torr 
CS2 produce essentially straight tracks typically 300-700 mm long [15] that 
are easily rejected from a dark matter analysis by range and energy cuts (see 
§5.3 and [13]). Such tracks provide an ideal means for developing track re- 
construction algorithms. Studies (see §4 and §6) of the alpha particle event 
populations found in DRIFT-II data also provide important information to 
allow elimination of such background events from present and future modules. 

The remainder of this paper will cover detector design (§2), data analysis (§3), 
a discussion of the different types of alpha particle events observed (§4), sim- 
ulations of alpha activity (§6), determination of the drift velocity and range- 
discrimination of different alpha particle tracks (§5). 



* Tel. +44-(0) 114-2223517; fax. +44-(0) 114-2728079 

Email address: t.lawson@sheffield.ac.uk (T.B.Lawson). 
^ The authors are aware of a new version of this code (SRIM2006), however this 
'includes no changes to the basic calculations of SRIM-2003'[15]. 
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2 DRIFT Detector Description 



The DRIFT project aims to make continuing improvements to the design 
of DRIFT detector modules. DRIFT-IIb is an upgraded version of DRIFT- 
Ila [12] in which the significant differences are a simplified arrangement of 
electronics and the replacement of copper fieldrings with stainless steel ver- 
sions. Briefly, each detector consists of a 1.5 x 1.5 x 1.5 m^ stainless steel 
vacuum vessel housing a 1 m^ dual low pressure TPC filled with CS2 vapour 
at 40 Torr. The two modules were operated at different mean drift fields (620 
V/cm in DRIFT-IIa, 581 V/cm in DRIFT-IIb). The two 0.5 m long drift 
regions share a central plane of 512 stainless steel wires of 20 /xm diameter 
spaced 2 mm apart. This arrangement forms a highly transparent drift cath- 
ode. Readout of charge deposited within the drift regions is via two 1 m^ 
multi wire proportional chambers (MWPCs), 50 cm either side of the central 
cathode. The MWPCs consist of grid planes of 512 stainless steel wires of 100 
fim diameter. The grid planes are spaced 1 cm either side of an anode plane 
comprising 512 wires of 20 /xm diameter. These are arranged orthogonal to 
the wires in the grids. The wire spacing in all the planes is 2 mm. Veto regions 
are instrumented around the edges of the wire planes. 

The central section of the anode and inner grid planes are each grouped down 
to eight outputs. As such, any track that extends a distance greater than that 
sampled by eight wires (> 16 mm) will reappear on the same output channnels 
at a later time in the event record (see [12] for details of DRIFT data format). 
This phenomenon was characteristic of all the alpha particle tracks observed. 
Since high energy alpha particles have a range greater than 300 mm in 40 
Torr CS2, the hit pattern on the eight outputs repeats many times. Figure 1 
depicts waveforms for an example alpha particle event. 



3 Analysis Procedure 

3. 1 Data Reduction and Parameterisation 

A dedicated analysis was developed to specifically select alpha particle like 
events from the data. Alpha particle tracks are sufficiently different from those 
of nuclear recoils that the analysis procedures described in Ref [13] would not 
have been appropriate. 

During data acquisition all electronics channel£i| had their signals digitised 
and recorded when any anode signal exceeded the hardware threshold (200 

^ A total of 36 channels comprising eight anodes, eight grids, a grid veto and an 
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Fig. 1. An example alpha particle track. Each trace is derived from signals produced 
on every eighth wire in the MWPC and represents the amplified charge signal plotted 
over time. The event record consists of 1 ms of pretrigger plus 4 ms of postrigger 
data (trigger at t = 0). Tracks that cross more than eight wires 'wrap around' and 
reappear on the original channel. The labels down the right hand side (in order from 
top) refer to: difference between grid and anode vetos; anode veto; grid veto; sum 
of anode channels; sum of grid channels; anode channels (x8); grid channels (x8). 

ADC points, where 2048 ADCs is equivalent to one Volt). Throughout this 
paper 'above threshold' means that the absolute value of the waveform ex- 
ceeded the absolute value of the threshold unless specified otherwise. Once 
acquired, three virtual waveforms for each MWPC (the grid sum, the anode 
sum and the veto difference — see Figure 1) were then constructed from exist- 
ing waveforms. All waveforms (real and virtual) then underwent a fast Fourier 
transform (FFT) to identify frequency components introduced by unwanted 
signal pickup (specifically, components around 50 kHz associated with the 
drift cathode power supply). The unwanted frequencies were notched out and 
the result inverse Fourier transformed. For each waveform, the baseline was 
calculated using the mean signal value in the pretrigger region between —900 
and — 100/is from the trigger time (a total of 1 ms is recorded prior to the 
trigger). The data were then smoothed using a Savitzky-Golay filter [14] using 



anode veto from each of the two MWPCs. 
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Table 1 

The main analysis parameters extracted from the raw data. 
polarity: +1 for positive-going signal-profiles, —1 for negative. 



height: Maximum absolute excursion of waveform while above 

threshold (in ADCs). 



voltage_weighted_time: 'Centre-of-gravity' in time, of a signal-profile. 



anode_area: 



Integrated area of all polarity 



1 anode signal-profiles. 



wires: 



Number polarity = — 1 anode signal-profiles with height 
> 45 ADCs, equivalent to the number of physical anode 
wires in the detector on which significant charge was de- 
posited. 



start-time: 



Start of the FWHM of the first polarity 
signal-profile in a track. 



1 anode 



end-time: 



End of the FWHM of the last polarity 
profile in a track. 



1 anode signal- 



a 4th-degree polynomial. This allowed the effective analysis threshold to be 
reduced to ~ 9 ADC points. 

The raw and smoothed data were reduced to a number of parameters that 
were stored in several tables within a relational database that could then by 
queried by subsequent analysis routines. The following is a brief description of 
the key parameters used. The names of all analysis parameters in this paper 
are refered to using italics. 



3.2 Analysis Parameters 

For each output channel, the voltage waveforms were scanned for pulses, i.e. 
sections which exceeded threshold. These sections are refered to as signal pro- 
files. A signal profile was deemed to start and end at the nearest points, either 
side of the peak, where the waveform crossed the baseline. The main param- 
eters extracted from the raw data are summarised in Table 1. 

Individual tracks were identified by grouping together negative polarity anode 
signal profiles in time order to form one 'track' entity. This algorithm relied 
on the assumption that tracks were well approximated by a straight line. 
This is a reasonable assumption for high energy alpha particle tracks. In a 
given event, tracks on each MWPC were recorded separately. Thus, a particle 
track straddling the central cathode and depositing charge in each side of 
the detector was counted as two separate tracks. These partial tracks were 
recombined at a later stage of the analysis. 
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The high charge-density of alpha particle tracks resulted in large pulses in 
the anode output waveforms. Thus a high threshold (45 ADC) was chosen to 
select these without contamination from sparks, nuclear recoils and gamma 
ray events. It was found, from a random selection of tracks, that reducing the 
pulse height threshold from 45 ADCs to 30 ADCs increased the mean x- and 
2;-components by less than 0.2%. Increasing the threshold to 50 ADCs resulted 
in the loss of 5 ± 2% of alpha particle tracks with the x- and z-components of 
the remaining tracks left unnaffected. 

The X component of a track (parallel to the wire plane, perpendicular to 
anodes), was calculated using the relation Ax = {wires — 1) x 2 mm. Grid 
signals were smaller than those on the anodes, requiring a different method 
for obtaining the track y component (parallel to the wire plane, perpendicular 
to grids): Ay = {At/dt — 1) x 2 mm, where At is the difference between the 
start-time and endJ,ime and dt is the mean time between successive negative 
polarity voltage_weighted_times in the grid waveforms occuring between the 
startMme and end_time. The threshold for accepting signal profiles on the 
grids was 9 ADCs. 



3.3 Track Definitions 

Tracks having wires > 8, distributed across all 8 anode channels (i.e. at least 
one signal profile per channel) were flagged in the database as candidate alpha 
particle tracks. A major part of the analysis involved the selection of alpha 
particle events that crossed the cathode from one side of the detector to the 
other without crossing a veto region. Such events were labelled 'Gold Plated 
Cathode Crossers' (GPCCs) and were selected with the following additional 
cuts: 

1. An alpha particle track, as defined above, appeared on both MWPCs, 

2. The event had no pulses on the veto difference line with height > 45 
ADCs. 

3. anode_area/ wires < 30 Vxyus. This cut eliminated occasional amplifier 
feedback — a problem that affected the highest amplitude events and was 
resolved in DRIFT-IIb with replacement preamplifiers. This cut alone 
retained > 87% of all candidate alpha tracks. 

4. anode_area > 30 Vx/is. This eliminated events such as sparks that de- 
posited charge on all readouts but had a combined charge less than that 
expected for alpha particle events. This cut alone retained greater than 
98% of candidate alpha tracks. 

5. start_time greater than -100 fis. This cut alone accepted greater than 96% 
of events and was introduced to remove events in which the pretrigger 
region contained pulses from the preceeding trigger. 
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6. end_time < 3000 /is. This provided a 1000 /is buffer between the track 
end-time and the end of the record, ensuring that only tracks with a mean 
voltage_weighted_time separation > 1000 fis, on a given anode, could evade 
this cut illegally (equivalent to < 0.1% of all candidate alpha particle 
tracks). 

In summary, none of the cuts mentioned above had a significant effect on the 
measured track range or derived drift velocity (see §5), since the dimensions 
of the tracks were not affected. However, cut 6 did constrain the orientations 
of the longest tracks in DRIFT-IIa, accepted as GPCCs, to be greater than 
~ 70° from the normal to the wire planes. 



4 Types of Alpha Particle Events 

4-1 Alpha Particle Event Parameters 

Alpha particles were expected in the detector due to decay of ^^^Rn (ri/2 = 3.8 
days, decay energy = 5.59 MeV) and subsequent alpha decays of ^^^Po and 
^^^Po. Radon progeny are usually produced as positively charged ions in CS2 
and so are attracted to the drift cathode where they plate out. Alpha decays 
from ^^°Po are ignored in this analysis since these occur below the long-lived 
^^°Pb (ri/2 = 22.2 years) in the ^ssu decay chain and subsequent event rates 
are insignificant. 

In many materials traces of ^^^Th are also present at comparable levels to 
those of 238-[j ^]^gj 'j]-^ig contamination gives rise to the much shorter-lived 
^^°Rn {Ea = 6.288 MeV, T1/2 = 55.6 s). The daughter nucleus from this decay, 
^^^Po, is also produced positively charged in CS2 and although its half life is 
relatively short (ri/2 = 145 ms) this is plenty of time for it to be drifted to the 
central cathode and be plated out on a cathode wire surface. There are two 
subsequent alpha decays (from ^^^Po and ^^^Po) plus two P decays with the 
decay sequence ending in stable ^°^Pb. 

A key issue in the decay sequence is the generation of recoiling progeny nuclei. 
Figures 2 and 3 illustrate the decay sequence from ^^^Rn to ^^°Pb and ^^°Rn 
to ^°^Pb, respectively, showing how interaction with a single cathode wire can 
allow low energy recoiling progeny to be observed without detection of the 
accompanying alpha particle, mimicking WIMP-nucleon interactions in the 
detector. 

From arguments and calculations presented in §5, the drift velocities in the 
DRIFT-IIa and DRIFT-IIb detectors are determined to be 59.3 ±0.2 (stat) ± 
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plate-out 



Pb-210 on wire-surface 




Fig. 2. Flow diagram of one decay cycle, starting with the alpha-decay of ^^^Rn and 
ending with the long-lived (tx/2 ~ 22 years) ^^'^Pb. The three shaded boxes denote 
events that produce a recoil with no accompanying alpha track observable. These 
events are termed radon progeny recoils (RPRs) and are an important background 
since they mimic the recoil signature expected for wimp-nucleon interactions. 

7.5 (sys) ms"""^ and 56 ± 1.7 (stat) ± 3 (sys) ms~^, respectively. A condition 
imposed by the analysis of alpha particle events required that ionization from 
an event completed its deposition in either MWPC within a time period of 
t = 3 ms. Combined with the drift velocity, this limits the z-component of any 
track recorded by an individual MWPC to be Zmax < 178 mm (DRIFT-IIa) 
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plate-out 



Pb-208 on wire-surface 




Fig. 3. Flow diagram, analagous to that depicted in Figure 2, starting with the 
alpha-decay of ^^'^Rn and ending with stable ^''^Pb. 

or < 168 mm (DRIFT-IIb). Modeling alpha particle tracks in 40 Torr CS2 
using SRIM-2003 [15] indicates typical ranges of 300-700 mm. These ranges 
are listed in Table 2. Since these ranges all exceed z^ax, only a subset of alpha 
particle tracks are accepted by the analysis. The fraction of events that are 
accepted is covered in §6. Tracks on a given MWPC that extend further in z 
than Zmax may still have been unvetoed but this fact could not be determined 
from the data. Such events that 'fell off' the end of the data record were also 
incomplete and so could not contribute to energy spectra. The effects of this 
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Table 2 

Energies and ranges of alpha particle tracks and resultant recoiling progeny expected 
in 40 Torr CS2 due to decay of ^^^Rn and ^^'^Rn and their respective progeny. Ranges 
and longitudinal straggling ('Long, strag.') are as simulated by SRIM2003[15]. Note 
that the polonium decays will typically occur on or near the surface of a cathode 
wire so, for example, an outgoing alpha particle will implant the resulting recoil 
within the wire and vice versa. 



Isotope 


Ea 


Range 


Long. 


Recoil 


-^recoil 


Range 


Long. 




(MeV) 


(mm) 


strag. 
(mm) 




(keV) 


(/im) 


strag. 
(/im) 




5.48948 


334 


13.4 


218po 


100.82 


577.91 


119.61 


218Po 


6.00235 


383 


15.3 


214pb 


112.33 


628.54 


129.62 


214p„ 


7.68682 


567 


23.3 


210pb 


146.64 


745.25 


149.99 


220Rn 


6.288 


413 


16.3 


216po 


116.5 


631.99 


129.08 


216po 


6.778 


464 


18.1 


212pb 


127.9 


682.27 


139.07 


212p„ 


8.785 


701 


30.0 


208p]2 


168.9 


818.14 


162.45 



and other cuts are discussed in 56. 



4-2 Alpha Particles from Decay of'^'^'^Rn, "^"^^Rn and uncharged Radon Progeny 

Because ^^^Rn, ^^°Rn and their uncharged progeny are the only effectively 
gaseous alpha sources under consideration this is the only case in which an 
alpha particle track could both start and end in the gas, under the assumption 
that dust contamination in the gas fill is negligible. Thus only these types of 
events could produce tracks that cross the cathode and satisfy the other GPCC 
criteria (see §3.3 and Figure 4). 

Radon and uncharged progeny events were characterised by (i) a clear, ex- 
tended track appearing simultaneously on both sides of the detector with an 
extent in the z direction less than z^ax] (ii) a lack of significant veto signals and 
(iii) simultaneous end times for tracks on each MWPC, corresponding with 
the point where the track crossed the cathode plane. In all cases, the changing 
mean peak height along the track follows the variation in dE/dx, producing 
a characteristic Bragg curve that reveals the direction of the original alpha 
particle. Figure 5 shows a representation of the Bragg curve for the track in 
Figure 4. The Bragg peak occurs in the left detector in this case indicating 
that the original decay occured on the right. The position of the central cath- 
ode is also shown, coinciding with a slight suppression of pulse area due to 
drift-field distortions around the cathode wires and charge division between 
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n \ 1 \ 1 r n \ \ \ \ r 

-1000 1000 2000 3000 4O00 -1O0O 1000 2000 3000 4000 

Time in microsec Time in microsec 



Fig. 4. A cathode-crossing alpha particle track due to radon decay (a GPCC event). 

the two sides of the detector. Figure 6 illustrates the positioning of a typical 
cathode crosser track within the detector, along with other classes of alpha 
particle event. 

4.3 Alpha Particles from Decay of^^^Po, "^^^Po, '^'^^Po and ^'^^Po 

Decay of ^^^Rn produces ^^^Po in the form of a positive ion approximately 
85% of the time in air [18]. Once produced, the uncharged ^^^Po fraction 
will be unaffected by the drift field and will thermally diffuse within the gas 
until it decays. The charged ^^^Po recoils, meanwhile, would be expected to 
quickly plate out onto the cathode wires [19]. One can expect subsequent alpha 
particles in the decay chain to have their origins on the cathode plane rather 
than in the gas. Some of the resulting tracks will stop in the gas without 
entering a veto region. This set of criteria could also be met by alpha particle 
tracks from radon decay that do not cross the cathode. An identical process 
occurs in the case of 220Rn progeny. Examples of this class of alpha particle 
track are shown in Figure 6 (C, D and E). 



5 Drift Velocity Measurements 

Knowledge of the anion drift velocity in the fiducial regions of the DRIFT-II 
detectors was crucial for accurate track reconstruction. To this end populations 
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2000 4000 

Drift time (/xs) 



6000 



Fig. 5. Plot of anode pulse area against drift time (proportional to extent in 
z-direction) for the alpha particle event shown in Figure 4. The time-sequence of the 
pulses on the right detector has been reversed and added to the end of the left-hand 
pulse sequence to allow the whole track to be plotted. The Bragg peak is clearly 
visible on the left. The vertical line marks the position of the central cathode. 

of alpha particles were observed and the drift velocity was measured via two 
techniques outlined in §5.2 and §5.3 below. 



5. 1 Field Simulations 



The drift velocity Vd of CS2 ions is related to the mean drift-field E hy = 
where is the ion mobility. Previous studies [6] using mixtures of CS2 
with Xe or Ar/CH4 demonstrated this linearity, however, other studies of ion 
diffusion in pure CS2 [7] revealed a slight deviation from the expected linear 
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Rn-222: 



Rn-220 




DRIFT 



DRIFT 



o 
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o 



Po-218 & Po-214: 



Po-216 & Po-212 



or 




500 



500 



Fig. 6. Illustration of different classes of alpha particle event. The 1 fiducial 
volume is shown bisected by the central cathode (dimensions in mm). Tracks that 
cross the cathode but remain inside the fiducial volume and have z components 
< Zmax (A) are counted as 'gold-plated cathode crossers' (GPCCs, see main text 
for definition) . These are used as a proportionate measure of the total radon event 
population (A, B and C, where B and C represent radon decay events that fail one 
or more GPCC selection criteria). Charged polonium decays (D and E) never cross 
the cathode, since, typically, they plate out on a cathode wire before decaying. Note 
that the higher range of polonium decays means that a smaller proportion of their 
total population would be accepted by the Zmax cut. 



relationship. Both these studies used considerably lower drift fields than those 
employed in this work and so it is helpful to extend our understanding to cover 
the regime in which DRIFT-II is currently operating. 

In order to fully understand the drift of charge into the MWPCs and sub- 
sequent signal development, detailed simulations of the electrostatics and 
avalanche process have been carried out. The results of this work will be 
covered in full in an upcoming paper [20]. In summary, the electric field was 
calculated in 3D using MAXWELL [21] finite element software in combina- 
tion with the GARFIELD [22] gaseous detector code. The purpose of these 
calculations was to allow direct comparisons of drift velocity between different 
modules operated with different mean drift fields. The simulations confirmed 
a uniform drift field in both detector modules. The drift field in DRIFT-IIa 
was found to be 620 V/cm. In DRIFT-IIb the corresponding value was 581 
V/cm. A simple linear dependance of drift velocity on drift field was assumed. 
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OUTER GR I D \ 
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-collimator 



Po-210 source 



Fig. 7. Setup for measuring drift velocity. A ^^''Po source was positioned as shown. 
5.304 MeV alpha particles from this source were collimated with a 250 fim thick 
polyester film with a nearly complete annulus cut out of it such that the path of the 
alpha particles subtended an angle of (45=b5)° to the normal to the wire planes. The 
mean residual range in the drift {z) direction from the inner grid was then 104.19 
mm, assuming the range of tracks was as predicted by SRIM2003. The mean drift 
velocity Wq was then deduced from the mean duration At over which the ionization 
from these events produced signals on the anodes. 



5.2 Collimated Alpha Particles 



Alpha particles from a thin-film ^^°Po source (alpha particle energy 5.304 
MeV) were directed into the fiducial volume of the DRIFT-IIb detector as 
shown in Figure 7. Alpha particles from the source with range R mm were 
collimated to an angle of (45 ±5)° by a 250 /im thick polyester mask attached 
to the back of the MWPC strongback (see Figure 7). On average the alpha 
particles stopped Az mm above the inner grid. By measuring the mean time 
At over which these alpha particle tracks deposited charge on the anode (grid 
to anode drift times are assumed to be small) the drift velocity was calculated. 



The distribution of alpha particle tracks within the x-y plane was investigated 
to ensure that reconstructed tracks matched the expected distribution. A con- 
cern was that the small size of the grid signals in early DRIFT-IIb operation 
would result in an underestimate of the track y-component Ay, since some 
signal-profiles would fall below the grid threshold and so not contribute to 
the measurement. This effect was found to reduce the apparent Ay measure- 
ment by approximately a factor of two. Since this implied a sytematic error 
a^y ~ Ay, events were selected in which the apparent ^/-component was less 
than 25 mm and the x-component was greater than 100 mm, to reduce the 
influence of this error on tracks reconstructed in the x-y plane. As such, the 
total number of events was reduced from 1410 to 241. 
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Table 3 

Variation in combined statistical and systematic errors on the value of R2xy as the 
maximum value of Ay is changed. In all cases, the cut Ax > 100 mm was also 
imposed. 



Ay cut 




aR2 (stat) 


(^R2 (sys) 


aji2 (combined) 


15 mm 


130.9 mm 


6.5 mm 


1.8 mm 


6.7 mm 


25 mm 


127.3 mm 


0.6 mm 


4.0 mm 


4.1 mm 


35 mm 


124.4 mm 


0.5 mm 


5.9 mm 


5.9 mm 



In [23] the accuracy of an earlier SRIM version (SRIM2000) was found to be 
within ~ 10%. In this work SRIM2003 [15] was used to predict alpha particle 
ranges in low pressure (40 Torr) CS2 for comparison with the data. Since the 
mean angle between the tracks and the normal to the wireplane was 45°, the 
mean projected range on this plane, R2^y, was assumed to be equal to Az. 
The collimator allowed a flat distribution of angles between (40 ±0.7)° and 
(50 ± 0.5)°, where the estimated systematic errors arose from ~ 1 mm errors 
in the collimator dimensions. The mean angle of tracks through the collimator 
was 6 = (45±0.18 (stat)±l (sys))°, where the statistical error is that expected 
for a random sample of 241 tracks passing through the collimator and the 
systematic error was derived from errors in the collimator dimensions of ±1 
mm. Based on these assumptions an independant confirmation of Az was 
extracted from the data and compared with values derived from SRIM2003. 
Figure 8 shows the distribution of R2xy for tracks selected to have Ax > 100 
mm and Ay < 25 mm, where: 

R2^y = ^{Axf + {Ayf (1) 



The mean of the values calculated using Equation 1 was determined as R2xy = 
127 ± 1 (stat) ± 4 (sys) mm where the systematic error is based on the 
assumptions that crAx(sys) = 2 mm and crA,;(sys) ~ Ay. The effect of changing 
the value of the Ay cut is shown in Table 3 where the overall error is found 
to be close to a local minimum when the cut value is set to 25 mm. As the 
Ay cut becomes more severe. Ay and thus o"Ay(sys) are reduced, however, the 
rapid drop in statistics causes an increase in (TAj/(stat). 

With reference to Figure 7 the residual z-component of a track is given by: 
Az = R2^ cot(e) = 127 ± 1 (stat) ± 6 (sys) mm, (2) 

where the error ctaz on Az is obtained via the usual law of error propagation 
taking into account the uncertainties in 6 and R2xy. 
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Fig. 8. Distribution of the 2-dimensional range R2xy calculated using Equation 1 for 
coUimated ^^''Po alpha particles in DRIFT-IIb. Only tracks where Ax > 100 mm 
and Ay < 25 mm were selected, thus limiting the effect of a/\y{sta,t) on the overall 
error. 

The mean alpha particle range R is given by: 

R = d sec{e) + R2^ csc{e) = 349 ± 1 (stat) ± 6 (sys) mm, (3) 



where d = 120 ih l(sys) mm is the distance from the source to the inner grid 
plane. The value for R should be compared to the SRIM2003 estimate of 
RsRiM = 317 ± 13 mm — suggesting a range correction-factor for typical alpha 
particle energies in 40 Torr CS2 of = R/Rsrim = 1-10 ± 0.05, consistent 
with the ~ 10% accuracy figure of [23], mentioned above. 

A total of 1410 alpha particle events were selected using the cuts described 
in § 3.3. This was reduced to 241 events by the Ax and Ay cuts mentioned 
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above. The mean drift time was At = 2252 ±23 (stat)±10 (sys)/is, where the 
statistical error was derived from the spread in the data and the systematic 
error is an estimate of the constraint on timing accuracy due to the shaping 
time of the amplifiers. Combining with the measured value of Az gives a drift 
velocity: 



The only significant operational difference in DRlFT-IIa, compared to DRIFT- 
Ilb was a slightly higher drift field (by a factor 1.067). Hence the drift velocity 
in DRlFT-lla was estimated to be Vd = (61 ± l(stat) ± 3(sys)) ms""*^. 

5.3 Background "^"^"^Rn Alpha Particles 

For this study, the intrinsic alpha particle background in DRlFT-lla was used 
rather than a deliberately placed source. Under the assumption that these 
events, and in particular GPCC events (§3.3), were randomly selected from an 
isotropically distributed population, it was concluded that all possible track 
directions were equally represented in a suitably large sample. However, an 
inspection of the x and y components of tracks revealed a systematic under- 
estimate of Ay, similar to the effect observed in DRlFT-llb data (§5.2), while 
the X and z components were found to have very similar distributions (see 
Figure 9). As in §5.2, this effect was found to reduce the apparent Ay mea- 
surement by a factor of two, implying a sytematic error a^y ~ Ay. To reduce 
the influence of this error on the full 3D track range a cut was applied that 
excluded tracks with Ay > 60 mm. An attempt could then be made to dis- 
criminate on range to select a single track species (for example, alpha particle 
tracks from ^^^Rn decay) from the data. As such, the total number of events 
was reduced from 3498 to 1031. 

Analysis of the GPCC events revealed that alpha particle tracks arising from 
the decay of four isotopes could be separately identified. Figure 10 shows a 
histogram of track range R for GPCC events from a DRlFT-lla dataset where 
tracks were selected to have Ay < 60 mm. The range was calculated according 
to: 



where the drift velocity = 59.2 ms~^ was chosen to produce the narrowest 
peaks and therefore the best separation in the distribution shown in Figure 10. 
The dominant peak in this case (within a region of interest (ROl) {345 mm < 
R < 395 mm}) was attributed to ^^^Rn. The mean value of the track range in 



Vd = Az/At = (57 ± l(stat) ± 3(sys)) ms 



-1 



(4) 




(5) 
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Fig. 9. Example scatter plot of x and z components of GPCC tracks taken from a 
DRIFT-IIa dataset. A drift velocity of 59.2 ms~^ was assumed to calculate Az. A 
circular arc of radius 371 mm is shown, representing the maximum range of 222Rn 
alpha particles. Note that the higher range of ^^^Po alpha particle tracks is indicated 
by an arc-like feature at larger radius (see main text for details). 

this region (containing 768 events) was R = 371.1 ± 0.3 (stat) ± 11 (sys) mm 
(the systematic error was calculated under the assumption of a ±1 ms~^ error 
in the assumed drift velocity). This range is larger than the SRIM-predicted 
track range listed in Table 2 by ~ 10%, consistent with the results in the 
previous section. 

The above range discrimination procedure was used as a cut to preferen- 
tially select 222f^j2 decay events from the GPCC population, since these were 
the most numerous and thus subject to the lowest level of contamination 
from other decay events. Specifically, tracks were selected with a dominant 
x-component by selecting only those tracks with Ay < 60 mm to reduce 
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Fig. 10. Histogram of the 3D range R for GPCC alpha events in DRIFT-IIa selected 
to have Ay < 60 mm. The peaks are attributed to alpha decays of ^^^Rn (at R ~ 371 
mm), ^^'^Po {R ~ 425 mm), ^^^Rn {R ~ 458 mm) and ^^^Po {R ~ 515 mm). Vertical 
lines indicate the mean peak positions. These ranges are 10% higher than those 
appearing in Table 2. A low-level continuum is also visible below R ~ 350 mm, 
attributed to alpha particles that traversed part of a cathode wire and subsequently 
escaped into the gas with degraded energy. 



the influence of underestimated Ay values. For each ^^^Rn event, the x and 
y components, Ax and Ay and deposition time At were measured with as- 
sumed systematic errors of cij; = 2 mm, cr^ ~ Ay, at = 10 /is. In addi- 
tion, tracks were then selected that fell within the range region of interest 
{345 mm < R < 395 mm}, calculated for the case where drift velocities were 
in the range {57 ms~^ < < 63 ms~^} and coinciding with the dominant 
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peak. The drift velocity vd for each track was then calculated according to: 



where the ^^^Rn track range R222 was treated as a free parameter and varied 
until the sampling error of the resulting distribution of Vn tended to a min- 
imum value. Using this procedure, the best estimate of the drift velocity in 
DRIFT-IIa, derived from '^'^^Rn alpha tracks, was 59.3 ± 0.2 (stat) ±7.5 (sys) 
ms~^ when R222 = 371 mm. The distribution of in DRIFT-IIa, determined 
with the best fit value of R222, is shown in Figure 11. For DRIFT-IIb, the cor- 
responding value was estimated to be vd = 55.6 ± 0.2 (stat) ±7 (sys) ms~^. 
These values are consistent with earlier results. 

There are several useful pieces of information provided by the plot in Figure 10; 

(1) The relative peak positions, in order of increasing R, match very closely 
with those expected for the alpha decays of ^^^Rn, ^^^Po, ^^°Rn and ^^^Po 
if a correction factor of ~ 10% is applied (see Table 2). 

(2) The fact that cathode-crossing alpha tracks are observed with a range 
consistent with ^^'^Po (and possibly ^^^Po) suggests that not all radon 
progeny plate out — some progeny remain suspended in the gas for long 
enough to decay and produce GPCC events. The most likely explanation 
for this is that some of the progeny are produced uncharged. 

(3) On this assumption, since each ^^^Rn decay results in exactly one ^^^Po 
decay, the ratio of the statistics in the two main peaks, adjusted for 
their relative observation efficiencies, determines what fraction fu of ^^^Po 
progeny are produced uncharged in 40 Torr CS2. 

(4) A similar calculation performed on the next two peaks can be used to 
determine the uncharged fraction of ^^^Po in 40 Torr CS2. 

The determination of the values of fu, as suggested by points (3) and (4) 
above, is discussed in §6. 

Table 4 summarises all the drift velocity measurements obtained in this sec- 
tion. It should be noted that the assumption aAy (sys) ~ Ay yields a con- 
servative estimate for the systematic uncertainty since the change in track 
range (and hence in vd) with the Ay cut is smaller than this uncertainty (see 



The track range measurements show a small deviation from SRIM2003 range 
estimates for alpha particles in 40 Torr CS2. A. Mangiarotti et al. [24] note 
inaccuracies in SRIM for heavy ions, mainly due to discrepancies in the elec- 
tronic contribution to the energy loss. Electronic stopping is by far the most 
dominant energy-loss mechanism for alpha particles in CS2. The authors of 
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Fig. 11. Distribution of drift-velocity measurements for alpha particle tracks from 
^^^Rn decay in DRIFT-IIa selected to have Ay < 60 mm. 



Table 4 

Results of drift velocity measurements. Values in [brackets] are derived values based 
on the drift field ratio between the two modules. 



Method 


VB [ms-i], DRIFT-IIa 


VD [ms^i], DRIFT-IIb 


Cohimated ^^°Po as 
^^^Rn Q background 


[61 ± 1 (stat) ± 3 (sys)] 
59.3 ± 0.2 (stat) ± 7.5 (sys) 


56 ±1.7 (stat) ±3 (sys) 
[55.6 ± 0.2 (stat) ± 7 (sys)] 



SRIM quote a mean deviation from experimental data of ~ 4.1% for the stop- 
ping power of He ions[15]. 
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6 Monte Carlo Simulation of Alpha Activity 



A Monte Carlo code was written to determine the total rate of radon decays 
Djin occuring in the detector given an observed rate g x D^n of GPCCs, 
where g is a. geometric factor dictated by the cuts outlined in §3.3 equal to the 
proportion of events that would be observed as valid GPCCs. From this, and 
a knowledge of the gas flow-rate F and total mass of gas rrio in the detector, 
the rate = k/g of radon emanation occuring within the vessel could be 
determined, where k is the rate of emanated radon atoms emitting an alpha 
particle that is then observed as a GPCC event. These quantities, for a given 
radon species, are related by the expression: 
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(7) 



where t is the elapsed time from the point at which radon starts to emanate 
into the vessel after initial evacuation and flushing with uncontaminated gas, 
Te = T1/2/ ln(2) is the characteristic decay time of the radon (where T1/2 is the 
half-life) and p is the characteristic flush time of the vessel, related to the flow 
rate by: 

p=^(l-l/e) (8) 



The typical flow rate for DRIFT-IIa operation was F = 0.118 kg/day and the 
corresponding flush-time was p = 3.15 days. At a pressure of 40 Torr mo = 588 
grammes for the 1.5 x 1.5 x 1.5 m^ vessel plus input pipework. 

The simulation modelled alpha particle tracks due to the decay of ^^^Rn, ^^^Po, 
^^°Rn or ^^^Po. The process was simplified by assuming all alpha particle tracks 
were straight and had ranges C = 1.1 times greater than those listed in Table 2. 
The range correction factor C was applied to reflect the slightly longer ranges 
observed in the data compared to SRIM2003 predictions (see §5.2 and §5.3). 

Tracks for each of the decays were simulated in separate runs to determine 
(i) the value of the geometric factor g for radon decay events for a given 
drift velocity and range correction factor C, (ii) the effect of varying the 
drift velocity, C and polonium uncharged fraction fu on the relative num- 
bers of radon and polonium GPCCs. The best matches to the observed ratios 
of event counts seen in Figure 10 were found when /^(^^^Po) = (22 ±2)% and 
fu^Po) = (10011,)%. 

All tracks were defined in terms of a start point and end-point where the end 
point was chosen to lie in a random direction but at a fixed distance, specified 
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Table 5 

Results from simulation of 10^ alpha decays. The entries are the numbers of events 
that satisfied the following criteria: I. Track confined within vessel; II. At least partly 
in fiducial volume (FV); III. Fully within FV; IV. In FV and also crossing central 
cathode; V. In FV and z < Zmax'-, VI. In FV, z < z^ax, crossing central cathode and 
> 8 anode hits on each side (GPCC); VII. GPCC with Ay < 120 mm (see §5.3 for 
the motivation for this constraint). 



vd (nis ^) 


Track category 


222Rn 


218pQ 


220Rn 


216po 


59.8 


I 


66652 


84327 


73093 


55665 


(DRIFT-IIa) 


II 


27847 


77736 


70679 


35852 




III 


15581 


48084 


27889 


11401 




IV 


3877 


1159 


673 


4577 




V 


8360 


17822 


8413 


4535 




VI 


1583 


291 


83 


740 




VII 


197 


33 


3 


52 


56.1 


I 


66751 


84309 


73152 


55643 


(DRIFT-IIb) 


II 


27900 


77714 


70725 


35931 




III 


15602 


48052 


27794 


11438 




IV 


3876 


1176 


651 


4616 




V 


7867 


16611 


7856 


4226 




VI 


1161 


261 


63 


645 




VII 


160 


25 


3 
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by the alpha particle range, from the start point. Radon start points were 
randomly chosen to lie anywhere within the vessel volume, as were uncharged 
polonium start points. Charged polonium start points were randomly picked 
from a one m^ plane coinciding with the central cathode. Table 5 lists the 
outputs from typical runs in which the uncharged fractions for ^ispo and 
2i6pQ vv-ere set at 28% and 100%, respectively. These fu values were chosen to 
produce the best match to observed ratios of radon to polonium progeny (see 
§5.3 and Figure 10). The value of g is given by the ratio of the number in line 
VI of the table to the total number of simulated alpha decays (10^). 

The values of the geometric factors in DRIFT-IIa for both ^^^Kn and ^^Of^^, 
determined from multiple runs of the code, were g = (13.64 ± 0.16 (stat) ± 
0.60 (sys)) X 10-3 and g = (9.63 ± 0.09 (stat) ± 0.57 (sys)) x lO'^, respec- 
tively, due to the cuts described in §3.3. The difference is due to the lower 
probability of the longer "^^^lin alpha particle tracks being fully contained 
within the effective fiducial volume. A similar difference is seen in DRIFT- 
IIb {g = (12.18 ± 0.11 (stat) ± 0.55 (sys)) x lO"^ for ^^^Rn, g = (8.62 ± 
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0.09 (stat) ±0.45 (sys)) x 10~^ for ^^°Rn) where the lower drift field reduced 
the fiducial volume in the z-direction. In all cases the systematic error takes 
into account uncertainties in the track ranges and drift velocities and the 
statistical error represents the spread in the results from ten identical runs. 

The rate of GPCC events observed in DRIFT-II background data was used 
to estimate the total emanation rate into the vessel from all components. 
Figure 12 shows data taken from a series of DRIFT-IIa runs spanning a period 
of about two and a half weeks in the middle of the 2005 data taking period. 
These particular runs were chosen because the vessel had been evacuated, 
flushed and refilled immediately before acquisition started. A fit to the function 
given in Equation 7 is also shown where the only free parameter was the 
reduced emanation rate k. A clear exponential rise is observed in accordance 
with the model described by Equation 7 with a best fit value of k = (13.6 ± 
1.5) X 10-3 Hz. 

Setting t = oo in Equation 7 allows one to determine the steady-state radon 
decay rate as Djin = 0.4 ± 0.04 Hz and thus the total steady-state emanation 
rate was = l.OiO.l Hz, assuming ^^^Rn formed the dominant contribution 
to the observed alpha activity. 



7 Conclusions 

An independant analysis, separate from that used to search for nuclear recoil 
events, has been developed to study phenomena associated with the alpha 
decays of radioactive contamination within the DRIFT-IIa and DRIFT-IIb 
directional dark matter detectors. 

The alpha particles produced as a result of the decay of ^^^Rn, ^^°Rn and 
their progeny (easily rejected in a dark matter analysis) have proved useful 
in determining a number of operational characteristics. The drift velocity was 
determined by two methods, one in each of the operational modules (see Ta- 
ble 4 in previous section). The first method, employing a directed ^^'^Po alpha 
particle source in DRIFT-IIb determined the drift velocity to be vd = 56 ± 1.7 
(stat) ±3 (sys) ms"^. The second method, using the intrinsic alpha particle 
background in DRIFT-IIa found vd = 59.3 ± 0.2 (stat) ±7.5 (sys) ms"^. 
Additionally, the range of alpha particle tracks was found to be an efficient 
discriminating factor between different alpha particle emitting species. This 
highlighted the possibility of some radon progeny being produced within the 
fiducial volume in an uncharged state and thus not being swept onto the cen- 
tral cathode by the drift field. Consequently, the decay of these uncharged 
progeny were able to produce alpha particle tracks that occasionally cross the 
cathode (GPCC events). Combining knowledge of the relative populations of 
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Fig. 12. Rate of GPCC events observed in DRIFT-IIa. Also shown is a fit to the 
data using the function given in Equation 7. 

radon GPCCs and those of their immediate progeny with results of a Monte 
Carlo simulation allowed an estimate of the fraction fu of polonium progeny 
produced uncharged. In the case of ^^^Po, the uncharged fraction was found 
to be fu = 22 ± 2% and for ^^^Fo, fu = 100l|j5%. This knowledge may prove 
useful in future more complex simulations of alpha emitting contamination in 
DRIFT-type detectors. 
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